VEHICLE BODY MOTION REALIZATION METHOD AND APPARATUS 



Technical Field 

The present invention relates to a vehicle body motion realization method 
and apparatus, and particularly to a vehicle body motion realization method and 
apparatus in which a vehicle realizes a predetermined vehicle body motion. 



Background Art 

Japanese Patent Application No. 2003-024177 discloses four-wheel 
individual steering and individual braking/driving control which maximizes vehicle 
body generating force in the process of realizing a desired direction in which vehicle 
body generating force acts, and a yaw moment. In this technique, a utilization factor 
of a road surface friction coefficient \x between a tire of each wheel and the road 
surface (hereinafter referred to as |i utilization factor) derives an optimum control rule 
which minimizes ^utilization factor on the assumption that the direction in which tire 
generating force acts is considered as a constant manipulated amount. 

Further, Japanese Patent Application Laid-Open (JP-A) No. 1 1-348753 
proposes control of braking force based on slip ratio target values of respective wheels 
for obtaining target yaw moment, vehicle longitudinal force and lateral force being 
calculated based on variations of yaw moment, vehicle longitudinal force and lateral 
force to small variations in the slip ratio between the respective wheels. 



Problems to be solved by the Invention 

In the apparatus disclosed in the aforementioned Japanese Patent 
Application No. 2003-024177, a yaw moment is assumed to be small compared with 
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necessary vehicle body generating force, and an optimum solution is derived based on 
a performance function corresponding to a ratio between a desired vehicle body 
generating force and a utilization factor of a road surface friction coefficient that is a 
friction coefficient between the road surface and a tire on the assumption that a 
direction in which tire generating force acts on each wheel is in the range of ±11/2 to 
the direction in which vehicle body generating force acts. 

However, at the time of inhibiting a spin occurring on, for example, a road 
surface of which friction coefficient \i is low, or the like, a required yaw moment 
becomes larger in relation to the vehicle body generating force, no optimum solution 
can be derived on the assumption that the direction in which tire generating force acts 
on each wheel is in the range of ±7i/2 to the direction in which vehicle body generating 
force acts, as described in the apparatus of the aforementioned patent application. 
Further, in the technique disclosed in the aforementioned JP-A No. 1 1-348753, a 
steepest descent calculation method using an inverse matrix of four rows and four 
columns in a braking force distribution logic for obtaining a slip ratio of each wheel is 
carried out, and a calculation load thereof is large. 

The present invention has been achieved in the aforementioned 
circumstances, and an object of the present invention is to provide a vehicle body 
motion realization method and apparatus, which allows realization of a predetermined 
vehicle body motion irrespective of respective magnitudes of a vehicle body 
generating force and a yaw moment. 

Disclosure of the Invention 

In order to achieve the aforementioned object, a first invention is a vehicle 
body motion realization method in which a direction in which tire generating force 
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acts is calculated for each wheel by means of a performance function including the 
direction in which the tire generating force acts, so as to realize desired yaw moment 
and vehicle body generating force for obtaining a predetermined vehicle body 
motion, and the vehicle body motion is realized by using the calculated direction in 
which the tire generating force acts, wherein the performance function is a 
performance function corresponding to a ratio between a sum of respective squares 
of the desired vehicle body generating force and yaw moment, and a utilization 
factor of a road surface frictional coefficient, which is a frictional coefficient 
between a road surface which is supposed to have a constant condition for each of 
wheels, and a tire. 

Namely, according to the aforementioned invention, a direction in which 
tire generating force acts is calculated for each wheel by means of a performance 
function including the direction in which the tire generating force acts, so as to 
realize desired yaw moment and vehicle body generating force for obtaining a 
predetermined vehicle body motion, and the vehicle body motion is realized by 
using the calculated direction in which the tire generating force acts. 

In this case, when it is assumed that the yaw moment is smaller compared 
with the vehicle body generating force, even if the direction in which the tire 
generating force acts is calculated for each wheel by means of a performance 
function corresponding to the ratio between the desired vehicle body generating 
force and the utilization factor of a road surface frictional coefficient that is a ' 
frictional coefficient between a road surface and a tire, so as to realize the desired 
yaw moment and vehicle body generating force, in a case in which the yaw moment 
is not smaller compared with the vehicle body generating force, the predetermined 
vehicle body motion cannot be properly realized even if the resultant force of tire 



3 



generating forces is used. 

Accordingly, in the present invention, as the performance function used in 
the present invention, a performance function corresponding to a ratio between the 
sum of respective squares of the desired vehicle body generating force and yaw 
moment, and the utilization factor of a road surface frictional coefficient, which is a 
frictional coefficient between a road surface which is supposed to have a constant 
condition, and a tire is used. 

In this manner, according to the present invention, a performance function 
including respective magnitudes of the vehicle body generating force and the yaw 
moment is used. Therefore, as long as the direction in which the tire generating force 
acts is calculated from the performance function and is used, the predetermined 
vehicle body motion can be properly realized regardless of a balance between the 
magnitudes of the desired vehicle body generating force and yaw moment. 

Further, the direction in which the tire generating force acts may be 
calculated for each wheel by carrying out linearization of restriction conditions of the 
desired yaw moment and vehicle body generating force and using the performance 
function. Thus, when the restriction conditions of the desired yaw moment and 
vehicle body generating force are linearized, calculating load can be reduced. 

Specifically, for example, when the aforementioned wheels are four in 
number with two front wheels and two rear wheels, the direction in which the tire 
generating force acts on each wheel is obtained by means of a pseudo-inverse matrix 
of two rows and four columns prepared from the two linearized restriction conditions 
of the desired yaw moment and vehicle body generating force, and from the 
performance function. 

A vehicle body motion realization method according to the second invention 
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comprises: calculating means that calculates a direction in which tire generating 
force acts for each wheel based on detecting means for detecting state amounts of a 
vehicle, and a performance function including a direction in which tire generating 
force acts, so as to realize desired yaw moment and vehicle body generating force for 
obtaining a predetermined vehicle body motion; and control means that controls a 
steering angle of each wheel and braking and driving forces of each wheel based on 
the direction in which the tire generating force acts, which is calculated by said 
calculating means for each wheel, a value of the performance function used in 
calculation of the direction in which the tire generating force acts, and a utilization 
factor of a road surface frictional coefficient, which is a frictional coefficient 
between a road surface and a tire calculated based on the desired yaw moment and 
vehicle body generating force. The vehicle body motion realization apparatus of the 
present invention has the same operation and effect as those of the invention as 
described above, and therefore, description thereof will be omitted. 

The calculating means of the vehicle body motion realization apparatus 
according to the present invention can be made to operate as below as is the case with 
the foregoing. Namely, the direction in which the tire generating force acts is 
calculated for each wheel by means of the performance function in which the desired 
yaw moment and vehicle body generating force are defined as restriction conditions. 
In this case, the direction in which the tire generating force acts may also be calculated 
for each wheel by carrying out linearization of the restriction states of the desired yaw 
moment and the vehicle body generating force, and also by using the performance 
function. For example, in a case in which the aforementioned wheels are four in 
number with two front wheels and two rear wheels, the direction in which the tire 
generating force acts may be calculated by means of a pseudo-inverse matrix of two 
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rows and four columns prepared from the two restriction conditions of the linearized 
desired yaw moment and vehicle body generating force and from the performance 
function. 

Brief Description of the Drawings 

Fig. 1 is a diagram showing a relationship between a critical lateral force 
and a yaw moment in a vehicle model with right-and-left two wheels. 

Fig. 2 is a diagram showing a relationship between a vehicle body 
generating force and a yaw moment when in a critical state. 

Fig. 3 is a diagram showing a vehicle model. 

Fig. 4 is a diagram showing a coordinate system corresponding to a vehicle 
body generating force. 

Fig. 5 is schematic structural diagram of a vehicle body motion realization 
apparatus according to an embodiment of the present invention. 

Fig. 6 is a flow chart showing a vehicle body motion realization program for 
realizing a vehicle body motion realization method. 

Figs. 7A to 7D are diagrams showing calculation results for directions in 
which generating forces act on respective wheels according to the embodiment of the 
present invention. 

Best Mode for Carrying out the Invention 

Referring now to the attached drawings, an embodiment of the present 
invention will be hereinafter described in detail. 

Initially, using a performance function corresponding to a ratio between a 
sum of squares of a vehicle body generating force and a yaw moment in an 
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embodiment of the present invention, and a utilization factor (^i utilization factor) of 
a road surface friction coefficient |li between a tire and a road surface, a principle of 
calculating the acting direction of a resultant of tire generating forces will be 
described. 

First, a description of a relationship between a vehicle body generating 
force and a yaw moment at the time of critical running will be given. 

There has been conventionally derived a control rule which maximizes a 
vehicle body generating force when a desired yaw moment and a direction in which 
a vehicle body generating force acts are given. However, there exists a trade-off 
relationship between the magnitude of vehicle body generating force and the yaw 
moment, and in a case in which a large yaw moment is obtained, the vehicle body 
generating force becomes smaller. To this end, when a yaw moment that exceeds the 
limit is required, it is necessary to consider as to which of the vehicle body 
generating force and the yaw moment should be prioritized and to adjust inputs 
(desired vehicle body generating force and yaw moment) to a control rule leading 
logic. Here, with the aim of utilizing for the aforementioned adjustment, a trade-off 
relationship between a vehicle body generating force and a yaw moment is derived in 
an approximate manner. 

Fig. 1 represents a vehicle motion model approximated with right-and-left 
two wheels. In this model, it is assumed that load movement of each of the wheels 
therebetween is ignored and that lateral forces of critical friction are outputted from 
both of the right-and-left two wheels. At this time, if braking and driving forces 
having a magnitude u are applied to each of these right and left wheels, a yaw 
moment Mz is given as below. 

Mz = 2/u (1) 
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wherein / represents a distance from the centroid of a vehicle to each wheel. Further, 
the resultant force F of lateral forces turns out to be the following expression. 



wherein \i represents the road surface |i, m represents the mass of a vehicle, and g 
represents a gravitational acceleration. From the relationship between them, the 
following expression (3) is given. 



The expression (3) is discussed on the assumption of no load movement or the like, 
but it represents the relationship between a vehicle body generating force, a yaw 
moment and the road surface \i. This expression (3) can also be considered as the 
extension to a four-wheel model by defining as an average distance between the 
centroid of a vehicle and a position of each wheel. Fig. 2 shows the relationship of 
the expression (3) wherein n= 0.78 and / = 1.5[m], and the relationship between the 
yaw moment and a maximized vehicle body generating force when a control rule 
derived based on a conventional logic (a direction in which a vehicle body 
generating force acts = 0, 7c/4, n/2, roll rigidity distribution = one by application of 
an optimum map) in the light of load movement, or even load nonlinear 
characteristics of a friction circle, is applied to a detailed model in light of load 
movement of four wheels, or even load nonlinear characteristics of a friction circle. 
The expression (3) represents the relationship derived on the assumption of no load 
movement, but it can be found that even in a state in which load movement or load 
nonlinear characteristics of a friction circle is still further considered, a relatively 
good approximation is obtained. 




(2) 




(3) 



It is assumed that as shown in Fig. 3, a direction 6 in which force applied to 
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a vehicle body acts as resultant force of tire generating forces acting on four wheels, 
and a magnitude Fimax (i=l: front left wheel, 2: front right wheel, 3: rear left wheel, 
and 4: rear right wheel) of a critical friction circle of each of the wheels are known, 
a direction in which tire generating force acts on each of the wheels (an angle qi 
formed by an acting direction of the resultant force and a single-wheel generating 
force) for uniformly maximizing the grip margin level of each wheel is obtained with 
desired yaw moment and vehicle body generating force being ensured. To this end, 
modeling of restriction conditions for ensuring desired yaw moment and vehicle 
body generating force is first carried out. When coordinate transformation is carried 
out in which a direction in which the resultant force of tire generating forces on four 
wheels acts is designated as x axis and a direction perpendicular to the x axis is 
designated as y axis, a position of each tire, (x, y) = (aj, bj), can be described as 
below (see Fig. 4): 
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T 

fl, = — cos^ - Lj^ sin 6 (4) 
2 



Tr 

02 = — ^cos 6 — Lj sin 6 (5) 



^3 = —COS0 ^L^sinO (6) 
2 



T 

a^= — ^^cos^ + L^sin^ (7) 
2 



T 

b^-^sin0 + Lf cos 6 (8) 
2 



Z>2 = — ^sin^ + L^cos^ (9) 



T 

b^^-^sin0.-L^cos0 (10) 



b. ^-^%\ne-L,cose (11) 
2 



Note that 0 is set to be 0 at the time of direct-advance acceleration and the 
counterclockwise direction is defined as the positive direction. Further, when a yaw 
moment to be produced around the centroid of a vehicle at a current time point is 
designated as Mzo, a vehicle body generating force is designated as Fq, and the \i 
utilization factor of each wheel is designated as y, the following restriction 
conditions would exist in the direction in which tire generating force acts on each of 
the wheels (to define the counterclockwise direction as the positive direction in 
relation to a direction in which the resultant force of generating forces acts, that is, 
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an angle q; formed by the acting direction of the resultant force of generating forces 
and a single-wheel generating force). 

T^iTmx sin 9i + }^2m^ sin + J^im^x sin ^3 + rF^,^ sin = 0 (12) 
7*'ia«« cosg, +}F2^^ cosgj +;*-3„, cos^j +r^4n»x cosg^ =Fo (13) 

- ai?f|m.x COS ?j - a2>**2nmx COS " «3;^^„bx COS - a^yF^^^ COS 

+ ^i;*^)»Bx sin 9, + b2j^2n^x sin ^2 + ^s^^^smax sin 93 + b^T^^^ sin ^4 = M.q (14) 

The aforementioned expression (12) can be described as follows: 

sin?4 =0 (15) 

Further, when the aforementioned expressions (13) and (14) are organized by 
deleting y from each of the expressions (13) and (14), the following expression is 
given. 

i>iFoF^^ sin + b^^F^F^^^ sin + b^F^F^^ sin + b^F^F^^^ sin q^ 

- (ajFo + M,o)Fi^ cos 9, - {a^F^ + M,o)/^j„^, cos q^^ 

- {a^F^ + M,o)^3max ^3 " (^4^0 + ^roKmax COS ^ 4 = 0 (16) 

The restriction conditions are resultingly described as in the expressions (15) and 
(16). 

When the magnitude of vehicle body generating force required is 
sufficiently larger than the yaw moment as in the conventional art, the direction in 
which tire generating force acts on each wheel is close to the direction in which 
vehicle body generating force acts, and it is considered that the relationship -7i/2<C 
qi<C7c/2 has been established. In this case, the problem of maximizing the grip 
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margin level, that is, the problem of minimizing the |i utilization factor can be 
replaced by the problem of obtaining the relationship -7c/2<Cqi<C7c/2 that maximizes 
the following expression (17) as a performance function (first performance function) 
based on the aforementioned expression (13) that is the restriction condition of 
vehicle body generating force. 
p 

= — = ^in«x COS 9j + ^2^, cos + COS + F^^, cos q^ (17) 
7 

Further, if the yaw moment Mzo is sufficiently larger than the vehicle body 
generating force Fq, the aforementioned problem can be replaced by a problem of 
obtaining the relationship aj — 7i/2 <C qi <C ai+7r/2 that maximizes the following 
expression (18) as a performance function (second performance function) based on 
the aforementioned expression (14) that is the restriction condition of the yaw 
moment. 

r M.0 M,o 

+ bEl^smq^ sin^2 + Mis-sin q^ + ^^sin q, 

Af^O ^rO ^zO 



*=1 VV^zO. 



where 



a;=-tan-'-^ (19) 



a. 



This condition represents that the direction in which tire generating force acts on 
each wheel generally coincides with the direction in which a desired yaw moment is 
produced. In this manner, a search region in the acting direction of tire generating 
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force to be obtained varies between a case in which the performance function is set 
based on the expression (13) that is the restriction condition of vehicle body 
generating force and a case in which the performance function is set based on the 
expression (14) that is the restriction condition of yaw moment. In the present 
embodiment, by taking this property into consideration, a search region adapted to a 
balance between the magnitudes of desired vehicle body generating force and yaw 
moment is set by combining these performance functions. 

The relationship between the vehicle body generating force and the yaw 
moment in a critical region is given by the aforementioned expression (3). From the 
foregoing, by achieving maximization of the following expression (20) serving as 
the performance function (third performance function), it is considered that the 
balance between the vehicle body generating force and the yaw moment can be 
evaluated in a balanced manner. 



jJjr^l^ (20) 

r 



The aforementioned expression (20) can be described below from the relationship 
between the expressions (13) and (14). 



= S^^^o^fma, COS9, +2]iV„««M,o(-a, COS +6,sing;) 

4 

= J max |f " ^0 - «f £0 )cos 9i + i*, M sin q, } (2 1) 

This problem can be solved by a technique of repeatedly using a pseudo-inverse 
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matrix. First, if the restriction condition of the expression (15) is subjected to 
Taylor's expansion, it is rewritten as below. 

A 

S ^^xnax {sin (fli — )cos q,o } = O 

That is. 



4 4 

2 9/0 ^/ = 2 i^^o COS - sin } (22) 

1-1 i«i 



Further, Taylor's expansion of the expression (16) is described below. 

4 

4 

= S + )^in.ax {cos q^^ - (i?- - )sin } 

That is, 

4 

Z ^0 COS ? ,0 + (^.- Fo + M ,0 )sin } • q, 

4 

= {^^f^ofero COS g,o - sin 9,0) + (^/^o + ^zoX^/o "n + cos 9^0)} (23) 

where qjo is qi in the previous step. If the expression (21) is approximated with a 
secondary Taylor's expansion, it can be described as below. 
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+ fe.M^ojsin 9,0 + cos 9,0(9. - gio)- (gi " 9io)^| 

= 2]^'«»xf-7f^^^o -«/J»*zo)-«>s9,o +*jM^oSin9«,}9,^ 
f-i L ^ 

+ f/ ^0 -'»f^io)- (9io cos 9,0 -sin 9,o)+6iWio(9,o sin g,o + cos 9,o))9i 



r 2^ 



1- 



9,-0 



cos 9,0 + 9,0 sin 9,0 



9io 



2 ^ 



sin 9,0 -9,-0 cos 9,-0 ! 



= ii ^f™, [- 1 ft '^0 - CO )• cos + fc,M,o sin 4,0 for -^if+Y, 



(24) 



where 



^ ^ (^^^0 -g.^co )- fao COS g,o - sin 9,0) + ^M,ofao sin + cos g.-p) 

j- COS ^^,0 + ^i^zo sin 



1- 



9f0 



COS qjo + q,Q sin 



-sm9io -gjocos^io 



f/^Fp -aiM.o)-((?,o cosg,o -sing,.o)+^,.M.o(gjosin^,-o +cosg,o)f 
2(/'Fo - fl,M,o j- cos + 2f>,.M^o sin 9,0 



Here, if the relationship 
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(/^Fo -ajM,o) cosg,o +b,.M,o sin^^o >0 (25) 
is given, according to transformation of variable expressed (26) given as below. 



A = V^m« i^'Po - «i^^o )• COS + fe,M,o sin q,^ \\q,- X, ) (26) 

the maximization problem of the expression (24) becomes a minimization problem 
of the expression given as below. 

Incidentally, the condition in the expression (25) turns out to be the sum of a 
calculation condition when a performance function is set as in the expression (17) 
based on the expression (13) that is a restriction condition of vehicle body 
generating force, 
cos > O 

that is, 

l^FQCOSqfQ>0 

and a calculation condition when a performance function is set as in the expression 
(18) based on the expression (14) that is a restriction condition of yaw moment, 
cos - sin ^ ^ 

that is, 

M,o (- fl,. cos g,-o + sin q^^) > 0 

and it can be expected that the aforementioned condition is established even when 
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the vehicle body generating force and the yaw moment have the same level of 
magnitude. At this time, the expressions (22) and (23) into which the expressions 
(15) and (16) are subjected to a primary Taylor's expansion, can be respectively 
described as below (see the following expressions (28) and (29)) from the relation of 
the following expression (27). 

g/ = y (/ , \^ i Pi+X, (27) 



II 2 \ ' 

/=) y (/ Fo - Of M^o r 9iQ + bi • ^xo sin 

4 

= S ^ina, kfO COS - Sin q,o " X. COS 9,0 } (28) 
^ V P'irn^ {^.^0 COS gfo + {oiFo + ,o)sin g,o} ^ 

4 

= S ^'»-» {^J'^o (910 COS - sin 9,0 - X, cos 9,0 ) 

+ (a,Fo + M,o X^ro sin + cos 9,0 - X, sin q,^)} (29) 

Accordingly, qi that satisfies the expressions (22) and (23) and maximizes the 
expression (24), that is, minimizes the utilization factor y of \i can be obtained by 
repeatedly calculating the following recurrence formula: 



9i 

^2 
^3 

.^4 





£,2 


En 


E^: 


+ 




.^21 




E2, 


E^. 







(30) 
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where 



• cos + ZjjMjo sin 
£ ^ V^<max • fe-P'o cos g,o + (g.Fo + M^o )sin q,^ } 

4 

G| - 2 ^/-ax k.o <^os ^ ,0 " sin - ^/ COS ^ ,0 } 



^2=2 ^im" {^.^0 COS - sin - X, COS g,o ) 

+ (a,.Fo + M,oX-?ro sin 9,0 + cos - X, sin g^.o )} 
In the foregoing, represents a pseudo-inverse matrix. If it is assumed that E is a 
full row rank, can be obtained as below. 
E^ = • (E • E^y^ 

The expression (30) uses the nature that calculation of a pseudo-inverse matrix 
satisfies the restriction condition and derives a solution Pj that minimizes the 
Euclidean norm given as below. 

4 

2 



1=1 



Next, a description will be given of a specific structure of the present 
embodiment. 

Fig. 5 shows a vehicle body motion realization apparatus for realizing the 
vehicle body motion realization method of the present embodiment. Front left and 
front right wheels of a vehicle 12 are designated as lOFL and lOFR, respectively, 
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and rear left and rear right wheels which are driving wheels of the vehicle are 
designated as lORL and 10 RR, respectively. The front left and right wheels lOFL 
and 1 OFR which serve as a coupled driving wheel and a steering wheel are steered by 
a power steering device 16 of rack and pinion type driven in response to 
manipulation of a steering wheel 14 by a vehicle operator via tie rods 18L and 18R. 

The braking force of each wheel is controlled by controlling a braking 
pressure of each of wheel cylinders 24FR, 24FL, 24RR and 24RL using a hydraulic 
circuit 22 of a braking device 20. Although not shown in the drawing, the hydraulic 
circuit 22 includes a reservoir, an oil pump, various valve devices and the like. The 
braking pressure of each wheel cylinder is, at a normal operation, controlled by a 
master cylinder 28 driven in accordance with a pressing operation of a brake pedal 
26 by a driver. Further, if necessary, the braking pressure is controlled by an electric 
control device 30 as will be described below in detail. 

These wheels lOFR to lORL are respectively provided with vehicle speed 
sensors 32FR, 32FL, 32RR and 32RL which each detect a wheel speed Vwi (i = fr 
(front right wheel), fl (front left wheel), rr (rear right wheel), rl (rear left wheel)). A 
steering column to which the steering wheel 14 is connected is provided with a 
steering angle sensor 34 that detects a steering angle (p. Further, the vehicle 12 is 
provided with a yaw rate sensor 36 that detects a raw rate r of the vehicle, a 
longitudinal acceleration sensor 38 that detects a longitudinal acceleration Gx, a 
lateral acceleration sensor 40 that detects a lateral acceleration Gy, and a vehicle 
speed sensor 42 detects a vehicle speed V. Incidentally, the steering angle sensor 34, 
yaw rate sensor 36 and lateral acceleration sensor 40 detect the steering angle, raw 
rate and lateral acceleration, respectively, of which positive direction is defined 
along the left turning direction of the vehicle. 
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As shown in the drawing, signals which indicate wheel speed Vwi detected 
by the wheel speed sensors 32FR to 32RL, a signal that indicates a steering angle cp 
detected by the steering angle sensor 34, a signal that indicates a yaw rate r detected 
by the yaw rate sensor 36, a signal that indicates a longitudinal acceleration Gx 
detected by the longitudinal acceleration sensor 38, a signal that indicates a lateral 
acceleration Gy detected by the lateral acceleration sensor 40, and a signal that 
indicates a vehicle speed V detected by the vehicle speed sensor 42 are inputted to 
the electric control device 30. Although not specifically shown in the drawing, the 
electric control device 30 comprises, for example, CPU, ROM, RAM, and an 
input/output port device, and includes a microcomputer having a general structure in 
which the aforementioned elements are connected together by a bilateral common 
bus. ROM stores therein a vehicle body motion realization program for realizing the 
vehicle body motion realization method. 

Next, the operation of the present embodiment will be described with 
reference to the flow chart of a vehicle body motion realization program shown in 
Fig. 6. 

In step 102, various state amounts of a vehicle body detected by the 
aforementioned various sensors are captured. 

In step 104, variable i for identifying each wheel is initialized. In step 106, 
the variable i is incremented. 

In step 108, the maximum frictional force Fimax for a wheel identified as the 
variable i is calculated. 

In step 1 10, it is determined whether the variable i equals the total number 
io of wheels (in the present embodiment, four) or not. If it is determined that the 
variable i does not equal the total number io, there still exists some wheel for which 
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the maxirQum frictional force Fimax has not yet been calculated. Therefore, the 
routine returns to step 106, and the aforementioned step (step 108) is executed. 

Conversely, if it is determined that the variable i equals the total number io, 
there exists no wheel for which the maximum frictional force Fimax has not been 
calculated. Then, in step 112, the function shown in the expression (30) is defined 
using the maximum frictional forces Fimax obtained by the aforementioned 
calculation, and other physical quantities. In step 114, the direction in which tire 
generating force acts on each wheel (qj) is calculated using the function (the 
expression (30)). In step 1 16, by means of the calculated directions in which tire 
generating forces act on respective wheels, and the like, the braking and driving 
forces and steering angle of each wheel are obtained. In step 118, each wheel is 
controlled based on the obtained braking and driving forces and steering angle of the 
wheel. 

As described above, according to the present embodiment, a performance 
function not dependent on respective magnitudes of the vehicle body generating 
force and the yaw moment, and therefore, as long as the direction in which the tire 
generating force acts, which is calculated from the performance function, is used, a 
predetermined vehicle body motion can be properly realized. 

Further, in the present embodiment, the direction in which the resultant 
force of tire generating forces is obtained using desired yaw moment and vehicle 
body generating force as the restriction conditions. Therefore, it is possible to 
realize the desired yaw moment and vehicle body generating force and obtain a 
predetermined vehicle body motion with a high degree of accuracy. 

Moreover, in the present embodiment, the restriction conditions of the 
desired yaw moment and vehicle body generating force are linearized to minimize 
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the utilization factor of longitudinal friction coefficient, thereby resulting in 
reduction of calculating load. 

In the present embodiment, the direction in which the tire generating force 
acts is calculated for each of the wheels so as to realize the desired yaw moment and 
vehicle body generating force. This also means that, under the restriction conditions 
for realization of the desired yaw moment and vehicle body generating force, the 
grip margin level of each wheel is uniformly maximized, that is, the |i utilization 
factor is minimized. 

In the embodiment as mentioned above, the pseudo-inverse matrix is 
simply calculated repeatedly, but the present invention is not limited to the same. 
Alternatively, a penalty function may also be used. 

A penalty function having a combination of an evaluation on a deviation 
between a performance function J(q) to be maximized and the restriction condition is 
defined in the following expression given as below. 

= 7^ PpF, (qli ^\Jn^ i^) (3 1 ) 

where JpyCq) is the left side of the expression (15), JpMzCq) is the left side of the 
expression (16), and p is a positive number. Here, the solution q of the expression 
(30), and the penalty function relevant to a search point qo in the previous step are 
calculated. If p(q)<^P(qo)j the recurrence formula (30) is recalculated with q=qo. 
Further, if p(q)>P(qo), the expression q = (q+qo)/2 is considered as the next 
candidate search point, and the penalty function relevant to q is calculated. 

In this case, if P(q) <C P(q)5 the recurrence formula of the expression (30) 
is calculated again with qo=q. Further, if P(q) > P(q), the calculation q = (q + Qo)/2 
is carried out repeatedly until the relation P(q) <C P(q) is satisfied, and the step of 
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constantly decrementing the penalty function is carried out, thereby achieving 
convergence of the recurrence formula. 

The convergence of the recurrence formula (qi = qjo) satisfies algebraic 
equations (22) and (23) into which the expressions (15) and (16) are linearized 
around qio, and means that q; which maximizes the expression (24) into which the 
expression (21 ) is approximated around qio becomes qi = qio. Assuming qj = qio in the 
expressions (22) and (23), the following expressions are given as below. 

sin ?jo + ^2max Sin + ^3n,ax sln q^^ + F^^^ sin = 0 (32) 

l^lWn^ sin + ^2^0^2max ^20 + ^3^0^3rnax Sin q^ + b^F^F^^^ SID q^^ 

- {a^Fo + M^o)^i««x cos q^o " (<»2^o -^M.q)F^^^ cos ^^o 

-(^3^0 ^A^zoKmax COS - (^4 ^^0 + zO )^4n«x ^OS =0 (33) 

Thus, it can be understood that qio satisfies the expressions (15) and (16) and 
becomes a local optimum solution that maximizes the evaluation function of the 
expression (21). 

Incidentally, the n utilization factor y is calculated from the value qi thus 
derived, based on the following expression (34). 

r^-. 1^^^ (34) 

The logic of the expression (30) can be solved within a domain in which the 
expression (25) is established, and therefore, an initial value of the recurrence 
formula also needs to be set within this domain. The expression (25) can be 
described as below. 
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{l^Fo - a,M .o)- cos +biM,o sin 9,0 

= V(/'Fo -flfM.o)' +(6iM,o)' -005(9,0 -«,)>0 (35) 

where 

J'"";" (36) 

Therefore, it is believed that the following expression (37) is an initial value of the 
expression (30). 



q,o=a,^t^ (37) 

This angle represents the direction in which the tire generating force acts on each of 
wheels when the desired vehicle body generating force and yaw moment are 
obtained by generating the tire generating force acting on each wheel in the direction 
in which the vehicle body generating force acts, which direction is expressed below, 

^Fo-'f (38) 

and also in the direction in which the yaw moment is produced, which direction is 
expressed below, 

^ifeo-— 

(39). 

Further, the {j, utilization ratio y having a value of 1 or more in the 
calculation based on the expression (34) represents that a target vehicle body 
generating force and a target yaw moment cannot be achieved under characteristics 
of tire generating force at the current time point. At this time, if the direction in 
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which tire generating force acts, which is calculated by the expression (30), is 
realized by the magnitude of a critical friction circle, the vehicle body generating 
force of Fo/y and the yaw moment of Mzq/y are obtained. This means that when the 
target values of vehicle body generating force and yaw moment cannot be achieved, 
the vehicle body generative force and the yaw moment are uniformly made smaller 
until respective achievable values thereof are both obtained. 

Figs. 7A to 7D are results of the directions in which generating forces act 
on respective wheels obtained based on the expression (30), as an example of 
calculation, when the direction in which the vehicle body generating force acts is 
expressed as 7t/2. In this case, the recurrence formula of the expression (30) is 
solved with I = 1.5 m. In the case shown in Fig. 7A, that is, when the target vehicle 
body generating force is larger compared with the target yaw moment (Fo=1000[N], 
Mzo=1000[Nm]) and the \i utilization ratio y of each wheel is 0.6585, and also in the 
case shown in Fig. 7B, that is, when the target yaw moment is larger compared with 
the target vehicle body generating force (Fo=1000[N], Mzo=l 0000[Nm]) and the \i 
utilization ratio y of each wheel is 0.4218, it can be understood that proper solutions 
can be obtained for both of the target values. Incidentally, in cases shown in Figs. 
7A and 7B, the following relationship is obtained. 
(pFo-ajMzo)cosqio + biMzosinqio > 0 

Further, even in cases shown in Figs. 7C and 7D, that is, even when the target vehicle 
body generating force and the target yaw moment are both set to be large values 
(Fo=10000[N], Mzo=±10000[Nm] (in Fig. 7C, Mzo=10000[Nm], and in Fig. 7D, M^o= 
— 10000[Nm]), it can be understood that appropriate solutions are obtained. 
Moreover, a domain in which the expression (31) is established, which domains 
corresponds to a search domain of the optimum solution, is indicated by dots in the 
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diagrams shown in Fig. 7. These domains vary properly in accordance with the 
magnitudes of vehicle body generating force and yaw moment both to be required. 
As a result, it can be understood that the logic for obtaining the direction in which 
generating force acts on each wheel does not need to be switched to another based on 
the magnitudes of the vehicle body generating force and the yaw moment, and the 
direction in which generative force acts on each wheel can be constantly calculated 
using the expression (30). 

Effects of the Invention 

As described above, the present invention has the effect that, since a 
performance function including respective magnitudes of a vehicle body generating 
force and a yaw moment is used, as long as a direction in which a tire generating 
force acts, calculated from the performance function, is used, a predetermined 
vehicle body motion can be properly realized regardless of a balance of magnitudes 
of desired vehicle body generating force and yaw moment. 
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